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bigskip 

A wave-function framework for the theory of the (e,e'N) reaction is 
presented in order to justify the use of coupled channel equations in the usual 
Feynman matrix element. The overall wave function containing the electron 
and nucleon coordinates is expanded in a basis set of eigenstates of the nuclear 
Hamiltonian, which contain both bound states as well as continuum states.. 
The latter have an ingoing nucleon with a variable momentum Q incident on 
the daughter nucleus as a target, with as many outgoing channels as desirable. 
The Dirac Eqs. for the electron part of the wave function acquire inhomo- 
geneous terms, and require the use of distorted electron Green's functions 
for their solutions. The condition that the asymptotic wave function contain 
only the appropriate momentum Qk for the outgoing nucleon, which corre- 
sponds to the electron momentum k through energy conservation, is achieved 
through the use of the steepest descent saddle point method, commonly used 
in three-body calculations. 
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I. INTRODUCTION. 



The commonly accepted starting point for writing down the (e,e'N) matrix element is a 
Feynman diagram. The Feynman diagram consists of two vertices: one in which a nucleon 
is detached from the target nucleus, and the other which represents the electron-nucleon 
interaction, which scatters the nucleon into the appropriate final state. In 1970 Gross and 
Lipperheide [0] wrote down such a matrix element, and Boffi and collaborators based 
their work 0] on this approach. Walecka and collaborators 0], as well as Donnelly and 
co-workers also start from this type of approach. Additional refinements were gradually 
introduced. The optical model, initially used to describe the nucleon-nucleus distortion 0, 
was later refined by using the Random Phase Approximation for the nuclear excitations , 
electron spin variables were introduced ||^, the relativistic nucleon-nucleus interaction was 
considered []^, and two-step processes (or coupled channel effects) were included |Q, 0. 
Many good reviews of the (e,e'N) formalism exist, the latest one being due to Kelly [|10|. 

The objective of the present paper is to formulate the (e,e' N) reaction in terms of a wave 
function expansion, rather than Feynman- diagram matrix elements This treatment permits 
one a) to understand the assumptions which underlie the coupled-channel treatment of the 
hadron-nucleus channels, which are commonly inserted into the Feynman matrix elements, 
b) to include the effect of the (e,e' N) channels into the treatment of dispersion corrections 
to electron-nucleus scattering, c) to provide a foundation for the treatment of the (e,e' 2N) 
processes. 

The wave function method described here is similar in nature to the expansions employed 
in three-body reaction theory [^, ||T2|, [|13|. In particular, the saddle-point method is used 
to obtain the asymptotic limit of the wave function for the two particles in the continuum 
(the electron and the emitted nucleon), thus making contact with the matrix element written 
conventionally for the Feynman diagram mentioned above. However, the present formulation 
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avoids the Faddeev treatment by neglecting the electro-magnetic interaction between the 
two particles in the continuum, since it is weak. In the extension to the (e,e' 2N) case, the 
interaction between the two nucleons in the continuum being strong, has to be included. 
This can be achieved by means by a method similar to the continuum discretized coupled 
channel (CDCC) treatment [|I^], of the breakup in deuteron-nucleus scattering, still 



without using the Faddeev approach. However, due to the complexity of the description of 
the two-nucleon correlations, both in the initial and final states, the extension to the (e,e' 
2N) case is left to a future paper. 

Our treatment is as follows. The overall wave function, which contains both the co- 
ordinates of the electron and the nucleons, is expanded on a basis set of nuclear states 
which contain as much of the correlations as needed. In the continuum, these nuclear states 
contain both ingoing and outgoing waves, even though the physical boundary conditions 
require that only outgoing nucleon waves exist asymptotically in the overall wave-function. 
This requirement is satisfied, as is explained below. The coefficients of the expansion of the 
overall wave function into this basis set are functions of the electron coordinates, which in 
turn obey inhomogeneous Dirac equations. The latter are solved by employing relativistic 
distorted Green's functions, as is described in Appendix A. This discussion is presented 
in Section 2 by means of a simple, one-nucleon example. The generalization to the many- 
nucleon case is done in section 3, and the coupled channel formulation is presented in section 
4. Complications due to the non-orthogonality between the bound states and one-nucleon- 
in-the-continuum-states do not arise at this stage because in the original expansion both 
the bound and the continuum nuclear states are eigenfunctions of the same hamiltonian. 
This argument is not invalidated by the subsequent replacement of many of the channels 
by effective complex interactions, such as by T-matrices or optical potentials. The non- 
orthogonality between the bound and the continuum optical model wave functions is thus 
an artifact of the optical model approximation, but does not invalidate the overall theory of 
the (e,e'N) reaction. 

The treatment presented here is schematic. The electron-nucleon interaction contains 



many intricate questions which have been addresses before but which are left out here. 
Also the details of the relativistic Dirac description of the electron, some of which are 
included in Appendix A, are left out in the text of this paper. Spin and orbital angular 
momentum information is also left out. In a subsequent paper by J. Kelly on the use of 
coupled equations in (e,e' p) reactions, some of the details are included exphcitly. The 
questions of principle addressed in this paper are: i) what set of basis states are suitable in 
the description of the wave functions involved, ii) how can one be sure that in the final state 
only the outgoing parts of these wave functions survive, and iii) what are the corresponding 
coefficients (the T-matrix elements) in the asymptotic waves. The realization that only 
outgoing waves survive for both particles is achieved naturally through the use of the steepest 
descent method in the evaluation of a saddle point in the integration space. Such approaches 
are common in the three-body literature, [1^ , [1^ but the extension to the coupled channel 
case requires additional thought, which is provided here. 



II. THE BASIC IDEAS. 

The basic ideas used in the present approach are illustrated in this section by means of 
a simplified example. In this model the "daughter" nucleus is described by an inert core, 
and the target nucleus consists of one valence nucleon bound to this core. The interaction 
between the electron with the core, and the nucleon with the core, is represented by two 
mean field potentials VE(re) and VN{rn), respectively. The interaction between the electron 

( 

and the hadron is denoted by Ven ■ It is commonly expressed in the form A^j^, where is 
the electron current and is the vector potential produced by the nuclear current. The 
two mean field potentials are real. 

The overall wave function is first expanded in a set of partial waves 

'^{fe,fn) = '^JMji^e,rn)yj,Mj{re,rn), 
^^^"^ J,Mj 

and each radial wave is expanded into a complete set of nuclear basis states (j) 
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i^{re,rn) = Xiire)(i)iirn) + / XkQire)(l)Q\rn)dQ. 
i=i 



(2.1) 



In the above equation spin and angular momentum numbers J, Mj should have been in- 
cluded, but, they already are left out. 

The discrete functions (pi{rn) are radial eigenstates with negative energies of the nucleon 
Hamiltonian 



2m dr"^ 



(r) = ef0,(r), t = l,2,...B. 



(2.2a) 



and the positive energy eigenstates form a continuum 



2m dr'^ 



2m 



<l)o\r) < g < oo. 



(2.2b) 



The functions 0q \r) have ingoing as well as outgoing waves. Asymptotically they are 



given by 
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rq \r) ~ e"*" sin(Qr + ^) = ^ {H^+\Qr) - H^-\Qr) e-^^") 



(2.3) 



The i/(±)'s are related to the outgoing and ingoing Hankel functions according to H^^'(z] 
izh^^Hz), and H^^Hz) = —izh^'^Hz) , where the functions h^^''^^ are defined in 10.1.1 of [116 



and ?7 is a nuclear phase shift determined by V]sr{r). The orthogonality properties of the 0's 
are as follows: 



<p'i\r)^^n\r)dr = ^S{Q-Q') 



(2.4a) 



P{r)dr = e-'^^^5iQ-Q') 



(2.4b) 



< 



/.{±)i 



>= ; < (j)i'\(f)i >= 6i^i'. 



:2.4c) 



In the above. 
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The orthogonality between the continuum states and the bound states follows from the fact 
that the potential V^v in Eqs. (2a) and (2b) is the same, hence non-orthogonality is no 
problem here. 

The functions x, coefficients of the expansion (2.1), are to be determined from their 
respective Dirac Eqs., as will now be described. The Hamiltonian of the system is 

H^H, + VE + Tr, + VM + v,n (2.5) 

where and are the relativistic and non-relativistic kinetic energy operators acting on 
the electron and nucleon, respectively. The equations for the electron functions x(^e) are 
obtained by multiplying the radial part of the overall Schrodinger Eq. [Ti — E)ip{re, r„) = 
on the left with either < 0i(r„)| or < 0^^(r„)|, making use of Eqs. (2.2) and (2.4), and 
integrating over r^j. One obtains 

{H, + Ve- Ef'^)xr{T,) = -B,{r,), (2.6a) 

{H, + Ve- 42)Xfcg(re) = -SQ(re), (2.6b) 
where the inhomogeneous driving terms B are given by 

Bi{re) =< 0>(fV(re,rn) > 

i?Q(^e)=<0^"^|^^(fV(re,r„)>. (2.7) 

If one inserts the expansions (2.1) for tlj(re,rn) into Eqs. (2.7), then Eqs. (2.6) become 
integro- differential equations for the functions Xi and Xk- However, with the exception of 
the ground state term Xg(re)0p(r„) (here i — g), which is present in the incident channel, 
all other x-terms are of order Ven, and hence the inhomogeneous terms B in Eqs. (2.6) are 
of second order in v^n, and they can be neglected as being dispersion correction terms. We 
follow this procedure here and replace ip{re, rn) by ipj 

i^{r,,r^) ^iJi = f^+\re)Urn), (2.8) 
6 



in the expressions (2.7) for the B's. Here f^'^> is a spinor which satisfies 

{He + VE-El^^)f^+\re)^0. (2.9) 

It has both incident and outgoing waves and describes the scattering of an electron from the 
core potential, 

Regardless wether the replacement (2.8) is made or not, the solutions of Eqs. (2.6) have, 
by construction, only outgoing waves asymptotically 

Xi ^ t}'^ exp(i^,) (2.10) 
XkQ ^Tjf^ exp(i^fcQ) 

where 

ei = kire + ^M2kir,) (2.11) 
&kQ = kqre + ln{2kQre). 

In the above expression the 7's are the Sommerfeld parameters, as explained in Appendix 
A, and the angular momentum phases l7r/2 have been left out. The factors T in Eqs. (2.10) 
are spinors which are obtained from the solution of inhomogeneous Dirac equations, as is 
discussed in Appendix A. 

The asymptotic expression for '0(re,r„) is obtained by inserting (2.10) into Eq. (2.1) 

B 

Hre, Tn) ~ ^i'^ ^MOi)Hrn) (2-12) 
1=1 

+ jf T^J exp(^,J - (ifW(Q^n) - H^--\Qrn) e"'^^) dQ. 

The electron and nucleon momenta in the expression above are related by energy conserva- 
tion, and that is what produces the correlation between the two wave functions. Defining 
the total energy as E, the electron energy as E^i^^ — (Hcky + {rrieC^Y, and the kinetic 
energy of the nucleon as either e^^^^^ or E^q^ — {hcQY /2m^ then the conservation of energy 
equations are 
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) + ef ) = ^{Hchy + {meC^Y + ef ^ = E (2.13a) 



+ "q"^ = \l(^^kQ? + Kc^)' + (^cg)V2m = (2.13b) 

from which it follows that dQ/dk = —{k/Q){mc^/ek)- Accordingly, the first term 
exp(i^fc^)i?*^+) in (2.12), which gives rise to expi(re/c + r„Q), has a stationary phase point 
at Qo, given by 

go/A:o = (r„/re)(mcVef^), (2.14a) 

and the corresponding integral has a non-vanishing value which can be evaluated by the 
saddle-point method [|l^, while exp{9kQ)H^~^ has no such point and does not contribute. 



Using Eqs. (2.13), and neglecting the rest energy of the electron in comparison to the total 
energy, one obtains 

go ^ irn/re)mcyhc (2.14b) 
hcko ^ E - (r„/re) ^7710^/2, 

where m is the nucleon mass. As described in Appendix A, the result is 

[tQ eMOk,) Y^H^^\Qrn)dQ (2.15) 

''"■v/2W7M7^ifexp(0fcJijW(Qorn). 



- 2 5/2 V-'"-^ /V'^W-fco 
' e 

The above result shows that asymptotically the wave function contains outgoing waves in 
both coordinates rg and r„ in spite of the fact that the basis set in r„ contained ingoing and 
outgoing waves, and further, that r„ and r^. are not independent of each other. Rather, the 
choice of their ratio determines the value of the electron energy, and hence the corresponding 
value of the nucleon energy which is to be observed. The amplitude of the emitted electron- 
nucleon pair, defined in Eq. (2.10), is determined by the T-matrix T^, which, according to 
Eq. (2.7) and Eq. (A. 9), is given by 
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+ / \ * 

Ut\r,)y [<pQo i^n)) v,ni^{r,r^)drjr^. (2.16) 

Here /^o is a Coulomb distorted electron spinor incident on the daughter nucleus with energy 
corresponding to the final state, and V'(rer„) is given by Eq. (2.8). In evaluating the integrals 
in Eq. (2.16) the two variables rg and r„ are independent of each other. This expression is 
of the same form as given by the use of Feynman diagrams. Once the asymptotic value of 
the wave function is known, the corresponding differential cross section can be evaluated as 
is explained in the literature [0, |ll9| . 

As long as ip{rf,rn) contains only bound state components for the nuclear wave function, 
such as in Eq. (2.8), the integral converges. If, however, also unbound components are 
included, then the integral diverges, a difficulty which has already been pointed out by Ma- 
haux and Weidenmiiller [^]. This is because now there are two particles in the continuum, 
whose interaction (through v^n in this case), extends over all regions of space. This difficulty 
can be seen in a three-dimensional example in which all continuum wave functions are re- 
placed by plane waves. For the case in which is given by Eq. (2.8), the three dimensional 
integral can be done without difficulty in momentum space and yields oc g^^0g(q + k^), 
where q = k^ — kj is the difference between the initial and final electron momenta, kj is the 
momentum of the nucleon in the final state, and (pg is the Fourier transform of the bound 
state wave function 0g/r„ of the nucleon. However if 0g/r„ is replaced by exp(zk^ ■ r„), then 
the corresponding result is Tk oc ^"^^(q + kf — k^) . This shows that the integrals in each 
partial wave component are ill behaved. Thus, if components of the wave function with two 
or more particles in the continuum is to be included in the overall wave function ip, (such as 
one electron and one or more nucleons) and if the interaction between these particles is to 
be included, then a different formulation is required. This is commonly done in three-body 
calculations by going to the coordinates of relative motion and center of mass motion of the 
two interacting particles, and will be the object of a future paper on the (e,e' 2N) reaction. 

The processes in which the incident electron first excites the nucleus by promoting the 
valence nucleon into an excited state and then knocks it out from there, can be included by 
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replacing il^irern) in Eq. (2.16) by the sum of terms on the right hand side of Eq. (2.1). The 
functions Xi(^e) with i ^ g would be obtained to first order in by the Green's function 
solution of Eq. (2.6a). These terms would include the effect of dispersion corrections on the 
(e,e'N) reaction, but only for the bound state excitations. 

In summary, this section exhibits the main ideas of how a continuum basis set for the 
emitted nucleon can be introduced, and the asymptotic electron wave function can be ob- 
tained by a saddle point method The expression for the (e,e'N) matrix element is identical 
to the usual Feynman diagram formulation, but transcends it if virtual nuclear excitation 
corrections are included 

III. GENERALIZATION TO THE MANY NUCLEON CASE, AND DISPERSION 

CORRECTIONS. 

The nuclear continuum wave function basis introduced in the previous section is un- 
realistic, in that the distorting potential due to the nuclear core, V^f, is real rather than 
complex, as would be required for an optical potential. This difficulty will now be removed 
by introducing internal degrees of freedom of the core nucleus, with the consequence that 
the inert potentials V^r and Ve are replaced by the microscopic interactions between the 
nucleons, which in turn calls for the coupled channel formalism. The orthogonality between 
the positive and negative energy nuclear states will still hold because these functions are 
eigenstates of the same Hamiltonian. The subsequent introduction of effective complex in- 
teractions, such as the optical model or effective coupling potentials, does not invalidate 
the orthogonality, however, inevitable errors are introduced when the non local effective 
interactions are replaced by approximate local ones. 

The target nucleus has ^4-1-1 nucleons with coordinates r„, r^^, ...r^, and the coordinates 
of the bound nucleons, ri,r2 ...r^, are collectively denoted by ^. Antisymmetrization will 
be ignored, the nucleon to be ejected by the incident electron will have coordinate r„, 
the electron coordinate is Tg, and, as before, spin degrees of freedom, as well as angular 
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momentum variables are ignored. 

The Hamiltonian describing the nucleons is 7i^+i(r„, the kinetic part for the relativis- 
tic Hamiltonian of the electron is again denoted by i?e(i"e) and hence the full Hamiltonian 
is 

H — He + Ve,A+l+'HA+l, (3.1) 

where 

Ve,A+l{re,rn,0 =Ve,A{^e, + V^^\^e, Yn) ■ (3.2a) 

and 

A 

^;e,A(re,0=$^^^^^Hre,r,), (3.2b) 

i=l 

is the electromagnetic interaction of the electron with the nucleons in the core. In the above, 
v^^\e,i) is the microscopic electromagnetic interaction between the electron and a nucleon 
which, as described in Section 2,. contains the Coulomb as well as current terms. At this 
stage no restrictive assumptions about the nature of the electromagnetic interaction are 
made. If, for example the interaction with the pion cloud (the so-called two-body currents) 
are included, then v^^\re, rj) does depend not only on the position and velocity of nucleon i 
but also on the coordinates and velocities of the other nucleons. Although this possibihty is 
implicit in the above expressions, it will not be written explicitly. Likewise, the interaction 
of the nucleon with the core nucleons is given by 

A 

VnA^n,0-^V^''\rn,ri), (3.2c) 

i=l 

where v^^\n,i) contains both nuclear and electromagnetic interactions of nucleon n with 
nucleon i. 

A set of eigenstates $(^+^) of the A -|- 1 nucleon Hamiltonian is now introduced. They 
are of two types. The first type corresponds to the bound states of the {A + 1) system 

Ha^M^'^'^M - 4^^^^*r^^(rn,0, (3.3a) 
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while the second set represents continuum states, which obey 

^^+i*^tr'"Vn, = '"Vn, o- (3.3b) 

The continuum states ^^q^^'^ are the complex conjugate of the states ^'Q^^'^\vn, For 
the latter one nucleon is incident with momentum Q on the daughter nucleus of A nucleons 
in a state of excitation iq. The outgoing channels of these states contain one, two or more 
nucleons or clusters of nucleons in the continuum. The nuclear energy €^q^^^ 

is composed of the kinetic energy of the free incident nucleon, the excitation energy ej^^ 
of the daughter nucleus in the state io, and the separation energy of the nucleon from 
the ground states of the target to that of the daughter nucleus. The recoil energy of the 
daughter nucleus is not written explicitly The bound states ^^f"^^^ and the continuum states 
^Q^it^ are eigenstates of the same Hamiltonian with different energy eigenvalues, hence 
they are mutually orthogonal: 

< $f |$(^^) >= (3.5a) 



< ^'r' >- (3.5b) 

In the above, the ket < ^q^j^^' '' | implies complex conjugation. 

The states $(^+^) are now used as a basis set into which the complete wave function is 
expanded. Indicating the angular momentum algebra only very schematically one obtains 
for the expansion 

*(re,r„,0 = f E 

^ J.Mj 

xyj,Mj{re,rn,i), 
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X/(re)$i^+^)(r„,0+ 



^0 



XfeQ,«o(^e)^Q,i(, 



(3.6) 



where J and Mj are angular momentum quantum numbers and 3^j,Afj(fe, f„,,^) are angular 
and spin functions. This expression is the many-body generalization of Eq. (2.1). 

The coupled equations for the electron wave functions will now be derived, while the 
coupled equations for the hadron functions will be discussed in section 4. For this 

purpose the notation will be simplified in what follows by replacing the subscripts / by a , 
the subscripts (Q,io) by (3, i.e.,. 

la'' p 

the combined subscript a and [3 by 7, i.e., 

E+/<i«-E. 

and the (— ) superscript on the continuum states $^^^^'~''(r„, ^) is suppressed. The replace- 
ment of a continuous variable by a discrete one is symbolic at this point, but it occurs 
naturally in numerical applications when an integral over momentum is discretized into a 
sum over momentum bins, as is done in the case of break-up in deuteron-nucleus scattering 



a El- 



Inserting the expansion (3.6) into the basic equation (Ti — E)'^ = 0, dropping the angular 
momentum subscripts, multiplying on the left with < ^^'^^^ \ , integrating over all nuclear 
coordinates, and making use of the orthogonality relations (3.5), one obtains the coupled 
equations for the functions x-y- For the discrete electron wave functions Xa the result is 

{H, - Ej^^ - Ui%) Xa = -J2 U^c^^X, (3.7a) 

where 

^S(^e) = < \Ve,A+A^^'^ > ■ (3.8a) 

The equations for the continuum wave functions Xi3 are 

(ife - + ^S) X/3 + E U';lXa + E f/S' X.' (3.7b) 



p' 
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where, 



f^S' = < Ve,A<-'^ > (3.8b) 



t^^i = < ^-'^ Ve,A+M^+^) > (3.8c) 
and the electron energies appropriate to each channel are given by 

E^-^) = E-e(^+i). (3.9) 

Contrary to what was done in Eq. (3.7a), in Eq. (3.7b) the interaction Ve,n between the 
two particles in the continuum was separated out and placed on the right hand side of Eq. 
(3.7b). The reason is that the matrix element of interaction between the two particles in the 
continuum will diverge, unless the wave functions in the matrix element are also distorted 
by that interaction. As noted before, the electron-nucleon interactions which occur in the 
matrix elements (3.8) are completely general, and include many-body currents. 

Equations (3.7) are the multi-nucleon generalizations of Eqs. (2.6). The sum over 7 in 
Eq. (3.7a), contains the electromagnetic coupling to the excited bound and continuum states 
of the target nucleus. With exception of the term which refers to the ground state of the 
target nucleus, {a = g) these terms represent dispersion corrections to elastic and inelastic 
electron scattering. Some of the couplings to the bound-state excitations have been included 
in the literature in the calculation of dispersion corrections to elastic electron scattering by 
means of explicit coupled channel equations [^, p2[, but explicit coupling to the continuum 



states, which contain the effect of the (e,e'N) reactions on the elastic or inelastic electron 
scattering, and which arise from the terms 7 = /3 in the sum in Eq. (3.7a), have up to now 
been ignored. Similarly, the second and third terms in the first line of Eq. (3.7b), with the 
exception of the term a = g , also represent dispersion corrections. 
The following approximations will now be made: 

a) The electromagnetic coupling between the electron and the nucleon in the continuum 
will be ignored. This means that the terms in the second line of Eq. (3.7b) will be set to 
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zero. Should it be desired to include this coupling, then the structure of the wave functions 
has to be modified, by making a transformation to coordinates of relative motion between 
the two particles in the continuum, and the coordinate of the center of mass of the two 
nucleons relative to the daughter nucleus. Such transformations are common in three-body 
calculations |T^, |T^, [pTSf . 



b) Dispersion correction terms will be ignored. This means that the terms on the right 
hand side of Eq. 3.7a) are replaced by 

As a result, Eq. (3.7b) together with assumption a) reduces to 

(if, - + f/^j) X, = -^S^,. (3.10) 

Since in Eq. (3.7b) one is dealing with continuum bins, it is not clear how valid it 
is to ignore the bins which are the closest neighbors to bin (3. This point needs further 
investigation. Once these terms are ignored, one obtains Eq. (3.10), which is of the form of 
Eqs. (2.6 b), with Bqire) given by 

BQ,^, = f/gx, = < $Si'"Ve,A+l$r'^ > fi^\re). (3.11) 

Hence the effect of the inhomogeneous term Bg i^ on can be carried out in the same 
manner as described there. In Eq. (3.11) the function Xg was replaced by fg^\ the solution 
of the Dirac equation for a electron elastically scattered by the target nucleus in its ground 
state, because, under the assumption of no dispersion corrections, Xg reduces to the latter. 
The functions ^^q^^^' ^ in Eq. (3.11) is quite general, but because of the choice of the (— ) 
boundary condition, in the limit when the electron and nucleon coordinates go to infinity, 
only the nucleon channel associated with the daughter nucleus in state io survives in the 
integral term in Eq. (3.6). Nevertheless, the presence of the other channels included in 
$0^^^' ■* affects the value of Bq^i^, as will be discussed further in the next section. 
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The assumption b) of ignoring dispersion corrections is not crucial to the present formu- 
lation. Their contribution to the (e, e'N) matrix element is obtained by replacing Eq. (3.11) 

by 

Bq,o = E < *Si!f'"Ve,A+ie^^^ > Xi+Hre), (3.12) 

a 

where to first order in the electromagnetic interaction the functions Xa are obtained from the 
solution of Eq. (3.7a), in which the sum over 7 is replaced by the single term U^^ fg^\fe)- 
In summary, the main point of the present section is to include the many-body degrees of 
freedom of the target nucleus by enlarging the meaning of the nuclear basis functions 
and ^^qX'~^ ■ The rigorous calculations of these functions is impossible at this point. 
However, the stage is now set for making practical approximations, by writing restricted 
sets of coupled equations for these functions, as is explained in the next section. 



IV. THE COUPLED EQUATIONS. 

The functions ^^'^'''^''(r^, ^) defined in the previous sections are general, containing zero, 
one, or more nucleons (or clusters) in the continuum. These functions are used in the 
expression (3.11) for the quantities Bq i^ , which in turn are required to obtain the asymptotic 
value of the electron and nucleon wave functions. In the present section an illustration is 
given of how to calculate the quantities -B^.^g when the functions ^ are restricted to 

a subset which contains not more than one nucleon in the continuum. The coupled channel 
treatment is analogous to the customary case of scattering of a nucleon from a target nucleus 
containing A nucleons. However, the present application differs from the conventional elastic 
scattering treatment in that the daughter nucleus can be in a highly excited state io , and 
hence the effective interactions, such as the optical model, will in principle be different from 
the case that io represents the ground sate. There are additional reasons why the use of an 
optical model is not appropriate under these conditions, as will now be discussed. 

The restricted set of functions ^if'^^^ will be denoted by 
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P$(A+1) 

where P is a projection operator which ehminates all components of which have more 

than one nucleon in the continuum. The complementary projection operator is Q , so that 
(P + Q)^^^^^ = The functions P^^^^^-* can be expanded in terms of the bound 

states ^['^^ of the daughter nucleus 



P<fi^+^nrn,0 =E-^'>.(r")'^J^^(^) (4.1) 



where the expansion coefficients (y9^^i(r„) obey the usual coupled equations driven by the 
nucleon-nucleon interaction and where the functions ^['^^ {^) are bound state eigenfunctions 
of Ha 

^A$!^^(0 = (4^^ + e.)$S^^(0, (4.2) 

Here ef is the energy of the state above the ground state, and is the separation energy of 
the nucleon from the ground state of the {A + 1) system to that of the A-body system. The 
total nuclear Hamiltonian is 

Ha+i = T^K) + Vn,AK. + ^a(0 (4.3) 

where T„ is the kinetic energy operator acting on the nucleon with coordinate r„ , and Vn,A 
was defined in Eq. (3.2c). 

The functions P$^/^^^'' obey the equation 

[PHa+iP - ^i"^) P$i^+^) = 0, (4.4) 



where |^| 



P7^^+iP = T„ + v^^JJ^ + PHaP, (4.5a) 
and where the effective hadron interaction is given by 

- Pv^^aP + PvuaQ (4^^'^ - Q'Ha+iQ + It) Qvn,AP, (4.5b) 
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and where 

E^"^ = e(^+^) - e'^^ - es. (4.5c) 

The interaction v^^^/'' differs from Vn^A in that it includes the effect of the channels with two 
or more nucleons in the continuum 

The coupled equations for the expansion coefficients are obtained by inserting the expan- 
sion (4.1) into Eq. (4.4), multiplying the equations on the left with < $j^'^'*(^)|, integrating 
over the coordinates ^ , and making use of the orthogonahty properties of the states ^^f'^ . 
The resulting, well known, coupled equations are 

dp 



2m dr. 



\ - 4? ) ^-^Arn) + E / f^ftVn, r':^^,Ar'n)K = (4.6a) 

n J it J 



where 



In the above the angular momentum terms such as Z(Z + l)/r^ have again been left out. 
The potential uf^ is non local and complex and depends on the state 7 in view of the 
presence of the Green's function in Eq. (4.5b). For negative energies (-E^^i ^ 0) ^Tv should 
become real and the states (^a^j decay exponentially at large distances. For positive energies 
> 0) the states (pp^i obey the (-) boundary condition, i.e., asymptotically for all states 
i these functions are ingoing, with the exception of state io, for which they are both ingoing 
and outgoing. For this case the energy Ej^^ can also be expressed as 

The quantities -Bq.jq, defined in Eq. (3.11), will now be expressed in terms of the functions 
(fi^^i. By inserting Eqs. (4.1) into Eq. (3.11), one obtains 

BQ,io - E < iVe,A+l) ^'G^'\rn,0 > /i+H^e), (4-7) 

i 

where the subscript P stands for {Q,io). This is the main result of the present section. It 
differs from the conventional DWBA expression 
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Bq,, =< /i-i(r„) (^e,A+i) *[f^'^(r.,0 > ft\re) 

in that transitions to the daughter nucleus in the final state io can take place by going first 
to a different state i through the electromagnetic coupling interaction, and from there to 
the final state via the strong hadron coupling effects. The probability for the latter step is 
given through the amplitude ipp^i obtained from the solution of the coupled equations (4.6a). 
One of these states i could be an important doorway state, whose presence would spread 
the (e,e'p) transition strength to neighboring states, as long as the latter are coupled to 
the doorway state [^. A formal justification for this point has apparently not been given 
previously. 

Examples of the use of Eq. (4.7) are as follows, a) The effect of coupling permits the 
inclusion of the charge exchange process, which consists of an (e, e'n) transition followed by 
a {n,p) process [§], 0. b) If the state io is chosen to be highly excited, which corresponds 
to a large missing energy, then the coupling to low-lying excitations such as hole states of 
valence nucleons, provides a quantitative description of how the nucleon, first emitted by 
knock-out from a valence orbit, excites the daughter nucleus on its way out of the nucleus, 
c) Coupling to excited states of the nucleon, such as a A excitation,. can describe the role of 
nucleon excitation in the (e,e'p) process. All of these effects would be missed in the DWBA- 
optical model formulation. 

The number of states included in the sums over i in Eqs. (4.6) and (4.7) can be excessively 
large. It is usual to replace most of the states in terms of a few "important" ones by means 
of expressions involving channel coupling Green's functions. As a result the states which are 
kept explicitly are then coupled to each other by effective interactions which are different 
from those of Eq. (4.5b) because they simulate the effect of a different channel space. They 
are non local, poorly known, and difficult to work with. In addition, the electromagnetic 
interaction fe,A+i in Eq- (4.7) will also have to be replaced by an effective interaction. 



^eff 



Ve,A+l (4-8) 



as is shown in appendix C, and as has already been discussed ||2^. The effective interaction 

19 



I'gjj resembles a hybrid nucleon-nucleon T— matrix. It is a mixture of electromagnetic and 
nuclear interactions, as shown by the occurrence in Eq. (4,8) of both the quantities v^^^J^ 
and Ve,A+i ■ It is also non-local, in view of the occurrence of the many-body channel Green's 
function ^^Thus, if only one state is used in Eq. (4.7), and that state is distorted by an 
optical potential, then the electromagnetic interaction has to be changed accordingly. 

In spite of the fact that many channels other than the one being measured can contribute 
to Bkg^ig in Eq. (4.7), the only channel which survives asymptotically in the overall wave 
function is the channel io, since this is the only channel which has both outgoing as well as 
ingoing waves. This follows if one inserts the asymptotic form of the electron wave function 

XkQ,io ~ exp(i^«^(re) 

and the nucleon wave functions 

into the integral term in Eq. (3.6), and uses energy conservation to correlate the electron 
and nucleon wave functions. The saddle-point method of evaluating the integrals over dQ 
then eliminates all terms but the ones which have outgoing waves in both the electron and 
nucleon coordinates, and only Hl^\rn)^if\C) survives. One regains the asymptotic result 
given by Eq. (2.15), with the only difference that the inhomogeneous term Bq^^, which 
determines the asymptotic electron wave function factor T^^ , is now given by the more 
general expression (4.7). 

V. SUMMARY AND CONCLUSIONS. 

A method has been described which both justifies and criticizes the common practice of 
replacing the final state wave-function in the usual Feynman matrix element for the (e,e' N) 
reaction by a coupled channel wave function. The method points out the approximations 
which are implicitly made when the final state wave-function is replaced by an optical-model 
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distorted wave, or when only a too restricted set of coupled equations are used. In these 
cases the electromagnetic electron-nucleon interaction has to be replaced by an effective 
interaction given by Eq. (4.8), which compensates for the approximations made in the wave 
function. This correction goes beyond the frequently applied Percy Factor correction, since 
the latter only makes up for the omitted non-localities in the optical model potential, and 
not, as Eq. (4.8) does, for the multi step transitions between channels, which precede or 
follow the electromagnetic transition. Difficulties which are due to the non-orthogonality 
between single particle bound states and optical model scattering states, do not occur in 
the coupled channel formulation, because the basic states of interest are eigenstates for 
positive and negative energies, respectively, of the same full hadron Hamiltonian. A further 
advantage of this approach is that it provides the theoretical framework for following the 
flow of flux from the initial channel (the ground state of the target nucleus) to the various 
segments of the inelastic channel space, into which the (e,e'N) strength is distributed. This 
information is needed to properly interpret the nuclear excitation aspect of the transparency 
information, obtained from (e,e' N) or (/i, /i' N) reactions [?]. 

The method developed here establishes a complete wave-function framework for describ- 
ing the electron-nucleus system, so that the coupling between channels 0, ^ in the theory 
of the (e,e' N) reaction can be justiflably included. The basis set of hadronic wave functions 
has ingoing as well as outgoing nucleon waves, similar to what is the case for wave func- 
tions used in the calculation of transition matrix elements By the use of the saddle-point 
method in the integral over the continuum basis functions in the asymptotic limit, the in- 
going part of the hadronic wave function cancels. The method is of the type commonly 
used in three-body calculations, but does not use the rigorous Faddeev formalism since the 
electromagnetic interaction between the two particles in the continuum is neglected. 

The formalism is applicable, with several important modiflcations, to the description 
of the emission of two nucleons from the target nucleus. The interactions between the two 
nucleons in the continuum is generally ignored , and in order to include them a procedure 



can be introduced which has strong similarities to the description of the effects of deuteron 
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breakup on deuteron- induced nuclear reactions fl^, |T5|. Experiments already point p8[ to 
the influence of the (e,e' 2N) process in the emission of one nucleon, and hence a theory for 
this process is timely. 

In summary, a derivation was given of the coupled channel wave-function description of 
the emission of a nucleon from a target nucleus by an incident electron. The derivation 
enables one to make contact with the theories of the optical model or of the coupled channel 
description of nucleon-nucleus reactions, and is generalizable to the case of the (e,e'2N) 
process. 
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Appendix A: 

Asymptotic solution of the inhomogeneous Dirac Eq. 

If the Dirac eq. describing the scattering on an electron from the electrostatic field of 
a nucleus, has also an inhomogeneous driving term, then the solution can be written in 
terms of a relativistic distorted Green's function acting on the inhomogeneous terms. The 
purpose of this Appendix is to present expressions for the asymptotic form of the result., in 
two different forms. 1. Involving the relativistic spinors which are the regular and irregular 
solutions of the distorted homogeneous Dirac equation, and 2. Involving the solutions of 
the Schrodinger-like solutions of the second order equations which are equivalent to the 
first order Dirac equation. This second result is of value when a comparison between the 
relativistic and non-relativistic treatments of the distortion of the nucleons in the (e,e'N) 
reaction is desired. 

Following the notation of Mott and Massey, p. 228, a Dirac spinor of a scattering 



wave function, when expanded into partial waves, has two sets of radial functions for the 
upper and lower components. They are denoted by Gi{r) and Fi{r) for one set of solutions, 
and (r) and F„/_i(r) for the other set. Each of these still multiplies a two-component 

spinor of the spherical harmonic angular functions, as indicated in Ref. If one defines 



new radial functions Q\{r) = Gx{r) x [kr) and J^\{r) = Fx{r) x {kr), with A = Z or — / — 1, 
then the coupled radial Dirac equations are 

a_g^{r)-{^--)T^{r)=0 (Ala) 
dr r 

a^J^,{r) + {^ + -)g4r)=0 (Alb) 
ar r 

where 

a±= [E ±mc^ ±V - e<^]/hc. (A.2) 

In the above, the quantum numbers k are the eigenvalues of the operator K = —{a ■ L + h), 
which assume positive or negative integer values. They are related to the angular momentum 
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quantum numbers / and the index A as follows. For k>0,A = — / — 1 and k = I. For k < 0, 
A = / and k = —I — 1. The vector spherical harmonics y^,^ which multiply the functions 
Q and JF are eigenf unctions of the operator K. The two upper and lower components of 
the Dirac spinors are given by r~^^fc(r)3^K./i and -r~^iJF_K^^3^_K^^, respectively [^]. In the 
definition of a±, E is the total energy of the electron, mc^ is its rest energy, V is the is the 
scalar part of the potential, which is added to the rest mass and which is present in the 
case of a nucleon scattering from the field of a nucleus, and e$ is the electrostatic Coulomb 
potential which is the vector part of the potential, added to the energy. The latter is negative 
(positive) if electrons (protons) interact with the positive charge distribution of the target 
nucleus. The asymptotic values of Q are 

^ sin(0z + ?7«) exp(i?7«) (A.3a) 

where 

(j)i = kr- ln/2 + 7ln(2fcr). (A.3b) 

The value of / in Eq.(A.3b) is related to k as described above, 7 is the Sommerfeld parameter 
Ze^/hv, and k is the wave number, related to the energy according to ^y {hcky + {mc^y = E. 
The inhomogeneous equations to be solved are 

a_P.(r) - (-^ - -)Q.(r) = fof )(r) (A.4a) 
ar r 

a+Q.(r) + (-^ + -)P.(r) = 6(^)(r) (A.4b) 
ar r 

Here V and Q replace the functions Q and JF in the homogeneous Dirac equation and fc'-^^ 
and b^^^ are the upper and lower components, respectively, of the inhomogeneous driving 
term. Asymptotically thes functions V and Q have only outgoing terms, the coefficients of 
which define the T-matrix elements Tp and Tq 

VKC^T'p^exp{i4>i) (A5a) 
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Q« ~ Tp^exp(i0,)- (A5a) 



Following the expressions given by Onley |3T[], for the relativistic Green's function, one 
obtains for the solution of Eq. (A4) the result 

V{r) = gR{r)A{r) + gi{r)B{r) (A.6a) 

and 

Q{r) = J^R{r)A{r) + Ti{r)B{r) (A.6b) 

where 

A{r)= {Gih^^^ +Tih^''^)dr' (A.7a) 

J r 

and 

B{r)= r {gnb^''^ +J'Rb^''^)dr'. (A.7b) 
Jo 

In the above the pair of functions {Qr-,Tr) are the regular solutions of Eqs. (A.l). with 
asymptotic behavior for Qr given by (A. 3a). Similarly, the functions (^/, jFj) are the irregular 
solutions, with Qi normalized such that the asymptotic value is given by 

Qi^ ^^exp(z0O- (A.8) 

The asymptotic behavior of the corresponding functions JF is obtained in terms of the Q 
's from Eq. (A. lb). The relativistic equivalent of the Wronskian between the regular and 
irregular solutions is 

^rQi—^iQr = 1 

whose value is a constant, which is equal to unity for the choice of the normalizations of 
^ijand Qi as described above. 

As a result of the above, the asymptotic value of the spinor {V,Q), defined in Eq. (A. 6), 
is of the form given by Eq. (A. 5), with the spinor 
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{E + m^)/{hck) 



\ 



(A.9a) 



V 



v 



I 



J 



and 



/ 

^0 



oo 



= S(oo) 



[QRh^''^ + Tnh^"^^) dr. 



(A.9b) 



This is the main result of the Appendix. 

The asymptotic expressions for the solution of the inhomogeneous Dirac equation will 
now be given in terms of the solution of the equivalent second order Schrodinger-like equa- 
tion. This procedure consists in transforming the two coupled first order differential Dirac 
equations for the radial partial wave functions (with or without the driving terms) into 
uncoupled second order differential equations, to which the usual wronskian techniques are 
then applied in order to extract the coefficient of the asymptotic outgoing wave. The " Dar- 
win factors" (akin to Perey factors) which are required in this transformation (in order to 
eliminate the first order differential terms) in the end cancel away in the term which re- 
sembles the non-relativistic expression, but remain for the additional terms which are also 
present. 

By differentiation and linear transformations, the first order coupled differential equa- 
tions (A.l) and (A. 5) are reduced to uncoupled differential equations of first and second 
order for T and Q, and similarly for V and Q. The coefficients in front of the Q" , Q' and Q 
are the same as the coefficients in front of P", V and V but in the equations for the latter 
the inhomogeneous terms and their derivatives also appear. Similarly for T and Q. The 
first order terms of these equations are then eliminated by introducing the Darwin factors 



F±(r) = 



(A.IO) 



and by performing the transformation 



(A.lla) 
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V{r) = F^{r)p{r) ; Q{r) = F^{r)q{r). 



(A.llb) 



one obtains second order differential equations for g and /, and similarly p and g, with 
unit coefficients in front of the second-order derivative terms. These equations for p and 
q of course contain the inhomogeneous terms, which are denoted hp and 6q, respectively, 
but otherwise they are identical to the equations for g and /. These transformations were 



already given by Darwin |2^, and are routinely used in the literature |^ in order to cast 
the Dirac equations into Schrodinger form. 

One obtains the asymptotic factor T-p, defined in Eq. (A. 5), by multiplying the equation 
for p" with g and the equation for g" with p, and integrating the difference over r from to 
oo. Because the coefficients in front of the g and p are the same, these terms cancel, and by 
making use of the asymptotic relations (A. 3) and (A. 5) one obtains 



I {gp" - Pg") dr = -kTp (F+(oo)) 
Jo 



-2 



(A.12a) 



The left hand side of the above equation is also equal to gbp dr. Combining these two 
results, one finally obtains, 

E + mc^ 



T 

■J- n 



+ 



pc 
E + mc^ 



g dr 



(A.13a) 



pc 



POO 






/ ^ 






Jo 




r J 



dr 



By proceeding in a similar way for and Q, one obtains 

/ {fq"~qf")dr = -ik-^J,{F_{^))- 
Jo E + mc^ 



-iTn 



(A. 12b) 



E + mc^ 

The left hand side of the above equation is also equal to fbg dr. Combining these two 
results, one finally obtains. 



Tg=i 



J^6(^) dr 



(A. 13b) 



POO 






/ ^ 


\{- 


--) 


Jo 







'/c 



dr. 



The first lines are the same as what one gets in the non- relativist ic treatment, however, 
the functions which enter into these expressions are the solutions ^ or JF of the Dirac 



27 



equation. If one used the non-relativistic simulations g or /of these functions, one would 
first have to multiply them by the normalized Darwin factors 



in order to transform them into the functions Q or JF. These factors are many times confused 
with the Perey damping factors, because they have approximately the same radial depen- 
dence [Q. However these two type of factors have a different physical origin and should 
not be confused with each other. The terms in the second lines in Eqs. (A. 13) are related 
to the spin-orbit potentials, since the term n/r is angular momentum dependent, and the 
factors a'^/a± are the coefficients of the spin-orbit potentials in the Schrodinger form of the 
radial Dirac equations. 
Appendix B: 

Asymptotic behavior via the method of steepest descent. 

In section 2 the wave function ^/'(re,r„) was represented in terms of an integral over a 
continuous set of basis functions 0q and the coefficients of the expansion, Xkqire), 

were obtained from the solution of an inhomogeneous Dirac equation with outgoing wave 
boundary conditions, exp(ifcQre). Asymptotically the function 4'Q\rn) has both ingoing 
and outgoing waves exp(±i(5r„), and hence the asymptotic form of ilj{re,rn) has terms of 
the form exp{ikQre ± iQvn)- The two continuum variables kg and Q are related by energy 
conservation, Eq. (2.5b). As Q increases, kg decreases, and hence ikgVe + iQvn has a 
saddle point at Qq where the integral over Q has a stationary phase The integral has a finite 
value which can be calculated by means of a steepest descent method around Qq. On the 
other hand, the integral with exp(zA;Qre — iQvn) has no stationary phase and it vanishes. 
The integral involving exp{ikgre + iQvn) can be obtained from the steepest descent method 
which gives [Arfken and Weber] the general result 



F±ir) = v^o^MTMoo) 




(B.l) 
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Here z — Q, 

s f{z) = ikqVe + iQ r„, (B.2) 

g{z) is a slowly varying function of a is a phase which depends on the direction of the 
path across the saddle point in the complex plane, and s is the variable which becomes 
asymptotically large. The saddle point is determined from the requirement df/dz = 0, 
which in our case gives 

^ = (B.3) 

where is the energy of the electron, as defined near Eq. (2.5), m is the mass of the nucleon, 
and X is the ratio 

x^-. (B.4) 
re 

Conservation of energy, Eq. (2.5b), was used in the derivation of Eq. (B.3). Prom the latter 
one can obtain Qq and ko separately. Neglecting the rest-mass energy of the electron relative 
to Ek one finds 

hcQo ~ mc^x (B.5a) 

and 

hcko ^E- mc^x^/2, (B.5b) 

where E is the combined energy of the electron and nucleon, as given in Eq. (2.5b). This 
result serves to fix the ratio of distances x. If the electron energy hckQ is determined from 
experiment, then the corresponding value of x is determined from Eq. (B.5b), and the 
corresponding value of Q from Eq. (B5.a). Inserting the value of 

d'^sf he 

into Eq. (B.l), one obtains 
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o • 2\ 1/2 



I{s) ^ g{Qo) — exp 

ncTp 



he 



(B.7) 



In the above, is assumed to be the independent variable, and r„ related to it via x. The 
exponent of Eq. (B.7) can also be written as exp [ireko + irnQo)] ■ 

Appendix C. Reduction of the number of channels in terms of effective in- 
teractions. 

The sums over i in Eqs. (4.6) and (4.7) contain a very large number of terms, which would 
be impossible to include in numerical calculations. Instead, it is more practical to eliminate 
most of the channels i in terms of a few leading channels by means of channel Greens 
functions and effective interactions p3[. . In this Appendix the resulting modifications to 
the (e, e'p) formulation will be briefly indicated. 

In order to simplify the notation, the potential vl^J will be replace by v, the matrix 
elements < ^i^^ {vl^A^l'f^^ > by < i\v\i' >, and the wave functions 0Q°j'~''(^n) will be replaced 
by (piivn). Two channels will be kept explicitly, and they will be denoted by a and b. It is as- 
sumed that the incident wave is in channel a. The kinetic energy operator — (/i^/2m) {(P/ dr"^) 
is denoted by T, and the channel energy Eg^i is denoted by Ei. The coupled equations (4.4b) 
then take the form 

(T — Ei)(j)i + < i\v\i' > (pii = — < i\v\a > 0^— < i\v\b > (pt, (C.la) 

i' 

(T+ < a\v\a > —Ea)(t)a + < a\v\i > (pii = — < a\v\b > 0^ (Cl.b) 

i' 

T+ < b\v\b > -Eb)(j)a + ^I^K' >(f)i' = - < b\v\a > 0„ (C.lc) 



Channels i can be eliminated formally p3[ by using the channel-coupled Greens function 



Gi,i'{rn,r'^), defined for example in Ref. By using these functions in Eq. (C.la) one 

obtains 

0i = ^ [< i'\v\a > (j)a+ < i'\v\b > 0b] , (C.2) 
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where integration over the second Green's function variable is implied. Upon inserting this 
result into Eqs. (C.la) and (C.lb), one obtains the "reduced" set of coupled equations 



(C.3a) 



where the effective potentials Uap are given by 



(C.Sb) 



) < i'\v\P . (C.4) 

Here a or P take on the values a and b. Because of the presence of the Gj^j/the potentials U 
are complex and non-local, and integration over the non-localities is implied in Eqs. (C.3). 
One can summarize the result of Eq. (C.4) by defining an effective nucleon-nucleon potential 



where the many body channel Green's function is given by 

If one makes use of the result 

(j)l<a\ + (l)l<b\ + J2(l>* <i\ = 

i 

4)1 [< a\+ < a\vg^] + (PI [< b\+ < b\vg^] 



(C.5a) 



(C.Sb) 



(C.6) 



one can rewrite the expression (4.6) for B in terms of that set of channel functions a and b 
which are kept explicitly as 



i:a~\rn) < 4T\^n) $r^(0 \viTf^'G^'\rn,0 > foiTe). (C.7) 



where 



^eff 



V 



(E) 
e,A+l 



(C. 
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This result shows that if only a few states are used in Eq. (4.6) for the calculation of the 
inhomogeneous term B, such as for example the single measured state |io > times the corre- 
sponding optical model wave function, then, in order to make up for the contribution from 
the other states an effective electromagnetic interaction has to be used . This interaction is 
a mixture of electromagnetic and nuclear interactions, as shown by the occurrence of vl^J^Q^ 
in Eq. (C.8). 
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